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Filter factors that determine an upper limit for the wave numbers for which refractom- 
eter measurements can be used to calculate the spectrum of refractivity are derived in this 
paper based on the assumption that refractometers measure a weighted average of the refrac- 
tive index in a volume of air surrounding the center of the refractometer cavity. Two models 
are used assuming the weighting function has spherical symmetry around the center and one 
model is used assuming the function has cylindrical symmetry. All models result in a simple 
mathematical form which should be easy to use in further theoretical developments. 

1. Introduction 

The wave number power spectrum of refractivity is determined for large wave numbers /■ 

by the correlation of the refractive index variations at points which are separated by small 
distances [I]. 1 Thus, refractometer measurements discriminate in two ways against the spec- 
trum for large wave numbers. First, the refractometer does not measure the variation of the 
refractive index at a point, but rather, ji weighted average of the variations in the volume of 
air in the vicinity of the center of the refractometer cavity, and second, the centers of two 
cavities may not be placed arbitrarily close because of obvious physical limitations. The first 
of these effects is evaluated in this paper assuming several different weighting functions and 
neglecting further practical limitations such as finite response time and noise characteristics of 
the instrument. The result shows that the measured spectrum is a filter factor times the 
theoretical spectrum. With appropriate approximations, the effect is also evaluated for the 
measured frequency spectrum of refractivity. 

The exact weighting function is determined by, among other things, the size and shape of 
the cavity and the size and shape of the openings which allow the passage of air through the 
cavity. The commonly used microwave refractometer cavity \2\ consists of a right circular 
cylinder with the ends partially opened to permit the flow of air and operated in the TE 0U mode 
at about 9 kMc Considering the weighting function for a closed cylindrical cavity operated 
in the TEqu mode, it appears that the weighting function for an open cavity with the same 
dimensions, and operated in the same mode, decreases approximately as a partial sine wave 
across a section of the cavity perpendicular to the axis and gaussianly along the axis from the 
center [3]. For mathematical simplicity, however, the most desirable weighting function is one 
that has spherical symmetry about the center of the cavity, as will be seen later in the analysis 
in this paper. (This would be a very good approximation, for example, for a refractometer 
using a cavity consisting of two concentric spheres, the outer sphere being a grid to allow for 
the passage of air.) Two of the weighting functions used in this paper have spherical sym- 
metry, and a third function is used which has cylindrical symmetry. The results for the three 
different forms differ significantly only in the region where data become unreliable because of 
other limitations. 

2. Theoretical Development 

It may be assumed that the refractometer measures a weighted spatial average of the 
variation An=n— n where w {) denotes an average value of the refractive index n. Let the 



I Figures in brackets indicate the literature references at the end of this paper. 
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weighted average of An at any time t be given by 



A v, ( I? 4\ *J 



An(R,t) = 



<PrA>t(r' ,t)f(r) \d 3 rAn(R +r, t)j(r) 



!: 



i\f(r) 



I 



d 3 rf(r) 



(1) 



where v is a volume centered at R, the center of the refractometer cavity, and/(r) is the weight- 
ing function for the cavity. Figure 1 shows two spherical cavities and the vectors associated 
with these. 




Figure I. Geometry for two spherical refractometer cavities. 
Using two refractometers with centers at R x and R 2 , it follows that 

( dV, | r?r a <An(r^OAn(ri,0>/(rO/(r 2 ) 
\ «' : V, f ^r 2 /(r,)/(r 2 ) 

J v. Jv 



<An(R l ,t)An(R 2 ,t)>-- 



where <C ^> denotes a time average or ensemble average. 
Following Wheelon [4] 2 , one obtains 



<AMr[,t)An&t)> = <An*>C(?[-/ i )=±J<PkS(k) exp [-it(r[-r' 2 )] 



(2) 



(3) 



where <^An 2 y> is the variance of the fluctuations of the refractive index. C(r[—r f 2 ) is the cor- 

-> -> -> 

relation of the fluctuations at two points determined by the vectors r[ and r' 2 , and S(k) is the 
wave number power spectrum of refractivity. 
Similarly, one may define 



<An(R u t)An(R 2 ,t)> = <An 2 >C w (P) 



(4) 



where C W (P) is the correlation function between the centers of the cavities of two refractom- 
eters, and the vector P is the vector between the two centers, P=R 1 — R 2 . 

2 The magnitude of the S(k) defined by Wheelon is about one-fourth as large as the one used here since he used Ae = 2nAn instead of An. 
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Combining (2), (3), and (4) one finds that 



<An(R 1 , t)An{R 2 , <)> <Ar, 2 >('JP)- 



[ ^ r 



dhr t <Aif>C(/ 1 -ii)f(rOf(r^ 



f dV, f (^/(^/(^[^JVfcexp [-tt-fa'-?,')]^*)] 
f d»r, f tPrJ^m 

'[ dhA dV 2 /(rJ/(r 2 )exp[-iMr 1 -r 2 )f 



ar'J, 



<Wfeexp[-' £fc-F]iSr(Jk) H 



f dV, [' 



cPr,f{ ri )f(r 2 ) 



(5) 



The factor in the braces in (5) is called the filter factor E and this determines the upper 

-> 

limit of k for which the refract ometer may be used to measure S(k). It should be noted that in 
the above formulation, the exact form of the structure of the atmospheric turbulence, whether 
isotropic or anisotropic, affects the weighting function but will have 4 no other effect on E. In 
practice, I* must be greater than the size of the cavity, but for theoretical purposes, one may 
assume that P can take on any value greater than or equal to zero. Then, taking the inverse 
Fourier transform of (5), one finds 

(6) 



ES(k)= cPP exp [ik-P]<An 2 >C w (P) 

which is the quantity obtained from measurements. ^ 

hi the following- we will evaluate E assuming the following forms for f(r): 



I for \r\<a, 



1°/W = - 



and 



LO for |r|>a 
2°/(r)=exp{ -r> 2 }, 



3°/(r) = - 



exp {— rljb 2 } iov ^r 2 -\-r 2 <a 
for V^+/1>a, 



where r X} r v , and r z are the components of r in the x, y, and z directions respectively. The 
first form, 1°, gives an upper bound for any closed spherical refractometer cavity of radius a, 
and the second form, 2°, gives an upper bound for such a cavity operated in the IEqh mode. 
In each of the three forms it is possible (but not necessary) to identify the parameter a (or b) 
with the dimension of the cavity. 

For the first case, 1°, if one has two refractometers which measure a small spherical volume 
of air of radius a and b, respectively 



r- sin XjrfXj d$i r\dr l exp [— ikr x cos \] 

(|«* 3 ) Jo Jn Jo 



x- 



a*) 



I-— sin X 2 ^X 2 dd 2 r\dr 2 exp [ikr 2 cos X 2 ] I 



(7) 
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The quantities in the braces in (7) are the same, except that the parameter (a) is changed 
and one is the complex conjugate of the other. Thus we may write 

E=g(k,a)g*(k,b) (8) 

where g*(k,b) is the complex conjugate of g(k,b). One then obtains by integrating 

g(k f a) = jj—r^ [sin ka—ka cos fa] (9) 



and thus, 



9 
E= „ \3/jL JA3 I sm ka—ka cos fa] [sin kb—kb cos kb\. (10) 



For the second case, 2°, one finds that 

E= < -3-372 sin \1d\1 dSi i\dvx exp [ — ikr\ cos XJ exp {— rl/a 2 } >■ 

X *< 13-^72 sm M^2 d0 2 r\dr 2 exp [ifa 2 cos X 2 ] exp {— r|/6 2 } V • (11) 

Again, one may set the two factors of the right hand side of (11) equal respectively to g(k,a) 
and g*(k,b). After performing the X and integrations, one obtains 

g(k,a) =-±= T s -^> exp { -r*/a>}dr (12) 

which yields, using relation No. 19, page 73, of reference^], 

, 7 4 a : %7r& f a 2 F\ f a 2 k 2 \ /1oN 

^ (M)= ^"4~* exp |-t-| =oxp \—r/' (13) 

Thus 

E=exp{-k 2 (a 2 +b 2 )/4}. 

Figure 2 shows the filter factors for the sharp spherical cutoff E s (ka) and for the gaussian 
cutoff E g (ka) versus fa, assuming that the refractometers used to determine E s (or E g ) are the 
same. 

For the case 3°, E\s given by 

E=g(a, b)g*(c, d) (14) 

where 

P dr x f dr y ("' *-, e- 7 ^' exp { -r//b*} 
g(a,b)= J -" J -' J -* (15) 

dr x \ drA r//'.exp{-r, 2 /6 2 } 

c/ — oo 1/ -a t/ — 1/ 



with L=^a 2 —r v 2 and #*(c, rf) is the complex conjugate of g(c, d). After integrating the de- 
nominator of (15) one obtains 

9(a, &) = ^ 2 T f *■. r dr, P dr z exp [-iirJb+rfo+rJti] exp { -r x 2 /b 2 }. (16) 

7T 5 <r0 J -co J -a J -L 
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Figure 2. Filter factors for gaussian cutoff and sharp cutoff. 




Figure 3. Geometry for locating the propagation vector k in 
relation to a cylindrical cavity with its axis parallel to the 
x-axis. 



The r x integral ion may he done using relation No. 23, page 121 in reference [5] and one 
obtains 



g(a,b) 



_ exp{-£/6 2 /4} f " , r 



ir<r 



.)>/ 



dr z exp {— i(r y k y +r 2 kz)}. 



Changing to polar coordinates, (17) becomes 

_exp {-k x 2 b 2 /4} C a 



<j(a,b)-- 



wa' 



ydy 

Jo Jo 



dd exp [—iy-yjk v 2 +k 2 sin 0], 



which may be integrated to yield 

g(a,b)-- 



t—Ir 2 h 2 ~\ 
~4 J 



a^'k y 2 +k z 2 
where J\(Z) is the Bessel function of order 1. Referring to figure 3, one finds that 

VV+&/=&Vl-sin 2 cos 2 jS 
and 

k x 2 =k 2 sin 2 8 cos L ' 



(17) 



(18) 



(19) 



(20a) 



(20b) 



where is the angle which k makes with the vertical, and (3 is the angle in the XZ plane with 
the x-axis assuming that the cavities are both alined along the x-axis. Combining (14), (19), 
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(20a), and (20b) it follows that 



E c - 



'ack 2 (\-sin 2 dcos 2 p) 



Jiiak-yjl—sm 2 B cos 2 p) J x {ck^\—%m 2 6 cos 2 p) 



X exp 



{ 



(sin 2 6 cos 2 P)k 2 (b 2 +d 2 ) 



} 



(21a) 



E c is plotted versus ka in figure 4, assuming two identical cavities, and also that a equals b in 
(21a), for the values P equal zero and several values of 6. 



One may also write 



Vl— sin 2 cos 2 /3=cos <£, 



where is the angle k forms with the z-axis. With this notation, (21a) becomes 
4 T , , T , , f k 2 (b 2 +d 2 ) sin 2 </>' 



^ 



r / z ^7^1 ^ / k 2 (b 2 +d 2 ) sin 2 </>\ 
A e/i (a& cos 0) e/i (d: cos 0) exp ^ —j- > • 



acF cos 2 </> 



(21b) 



In figure 4, the curves for equal to 0°, 45°, and 90° correspond to <j> equal to 90°, 45° and 0°, 
respectively. 

If one assumes isotropic turbulence and then carries out the angular integrations in (5) 

using the filter factor (21a), a new filter factor is obtained which, as in the two previous 

-> 

forms, is dependent only on the magnitude but not on the direction of k. 
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Figure 4. Filter factor for a cylindrical refractometer cavity 
with its axis in the horizontal direction. 

-> 
9 is the angle a vector k, in the xy plane makes with the vertical axis, and the 
-> 
angle 4> is the angle a vector k makes with the axis of the cylinder. 
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It may be noted here that zeros of the functions shown in figure 2 and figure 4 can actually 
occur, the only requirement being that the weighting function for the refractometer decrease 
slowly over a small volume surrounding the center of the cavity. 

The frequency spectrum of refractivity is given by [6] 

W n (p)= | dr cos vr<An(r u t)An(r u t+r)>. (22) 

Many measurements of this type, taken with a stationary refractometer, are available [2] 
Assuming (somewhat arbitrarily) that the structure of the turbulence is frozen in blobs which 
are carried past the refractometer by the wind with the const ant velocity ~ji, one may write 



-> -> 



<A<r„0An(r lf / + r)> <An(r 1 ,0An> 1 + M r,0>, (23) 

and then, using (1), it follows that 

<Aw(r 1 ,0An(r 1 , ^+r)>=<A^(r 1 ,OAn(r 1 +Mr,0>. (24) 

It should be noted that in the above expressions ((23) and (24)) no allowance is made for the 
"self motion" [4] of the air, and this may be of considerable importance. Applying (24) and 
(5) in (22), one obtains 

W na (v) = [ dr cos vt J , ( d?k exp [-ik • fXT] S(k) E(ka) , (25) 

J OTt J p 

where W na {v) is the spectrum which is actually measured and v=2wf where / is a frequency. 
For the case of isotropic turbulence (25) yields 

W na {v)=-^-\ dkkS(k)E(ka), (26) 

47T/XJ „/„ 

using either 1° or 2° for the filter factor. 

Substituting v=kp in (26) and differentiating both sides with respect to k } one obtains 

Norton [7] and Wheelon [4] each obtained the result 



(29) 



which, when combined with (27) yields 

dW n (kn) = 1 dWjky) 
d(kfx) E(ka) d(kn) 

which relates the slope of the actual spectrum to the slope of the measured spectrum. 

3. Discussion of Results 

Because no data are available for calculating the wave number spectrum, the interpreta- 
tion of results will be confined to t he frequency spectrum. It is well to bear in mind, however, 
tin 4 assumptions made (eqs (23) and (24)) to obtain these results. 

Referring to figure 2 and figure 4 4 it is seen that E(ka)^l for &«<0.5, so that in this range 
the filter factor has very little effect. 1 f one assumes that a=2.5 centimeters, which corresponds 
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to a micro wave refractometer cavity [2], it follows that k<^0.2 per centimeter. For the con- 
stant wind velocity n of 5 meters per second (approximately 11 mph) this yields /< 100/2 -k 
cycles per second as the range in which the frequency spectrum is unaffected by the filter factor. 
The available data [8] ^ives a spectrum accurate to frequencies of at most 10 cycles per second, 
and the filter factor does not alter the spectrum in this range. 

If the filter factor does go to zero, as in figure 2 and figure 4, the first zero determines a 
definite upper limit for k (or /) for which the refractometer determined spectrum is valid. 
However, refractometer measurements are generally inaccurate because of noise, etc., before 
this limit is realized. 

The need for more data to determine the weighting function and the two spectra is ap- 
parent. Even if such data were available, it is felt that one of the three forms derived here 
would give a satisfactory approximation. 
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